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We study the electron spin decoherence of encapsulated atomic hydrogen in octasilsesquioxane 
cages induced by the 1 H and 29 Si nuclear spin bath. By applying the Carr-Purcell-Meiboom-Gill 
(CPMG) pulse sequence we significantly suppress the low-frequency noise due to nuclear spin flip- 
flops up to the point where a maximum T2 = 56 [is is observed. Moreover, dynamical decoupling 
with the CPMG sequence reveals the existence of two sources of high-frequency noise: first, a 
fluctuating magnetic field with the proton Larmor frequency, equivalent to classical magnetic field 
noise imposed by the X H nuclear spins of the cage organic substituents, and second, decoherence 
due to entanglement between the electron and the inner Si nuclear spin of the cage. 
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Hybrid electron-nuclear spin systems and in particular 
paramagnetic atoms trapped in molecular cages like en- 
dohedral fullerenes (e.g. N@Cgo or P@Cgo) are promis- 
ing components of spin-based quantum computing be- 
cause they have long spin relaxation times and can be 
precisely placed into large arrays by chemical engineer- 
ing [l|, |2| . Examples of recent progress include storage 
of quantum states over 100 ms [3j, as well as design 
or preparation of new molecular architectures for scal- 
able quantum computing ka-|6J]. Atomic hydrogen can 
be even more attractive due to its simpler electronic Is 
state and the exceptionally large hyperfme coupling of 
1420.406 MHz. Systems with large hyperfme couplings 
can under certain conditions benefit from fast manipu- 
lation of nuclear spins as was recently demonstrated in 
bismuth-doped silicon p|. 

Whilst Ceo cannot stably host atomic hydrogen, it has 
been found that polyhedral octa-silsesquioxanes (POSS) 
are ideal cages for this purpose [8J|. Crucial properties 
for quantum computing like the spin-lattice T\ and spin- 
spin T2 relaxation times depend strongly on the type of 
the peripheral organic substituents [9(. Recently [lfjj, 
we showed that the room-temperature phase memory 
time Tm=13.9 /is for the species with R=OSiMe2H is 
the longest observed so far for this kind of cages. 

A key step in implementing a physical system as qubit 
is the identification of all possible sources of decoher- 
ence. Understanding the origin of magnetic interactions 
between the central spin and its environment can help to 
improve the performance of the qubit by prolonging the 
electron spin coherence using optimal dynamical decou- 
pling [111]. In addition, dynamical decouplin g te chniques 
allow for measuring noise spectral densities [12|, [l3[ and 
can thus map the electron spin environment. In this Let- 
ter we apply the CPMG pulse sequence on atomic hy- 
drogen encapsulated in POSS cages (H@POSS) in order 
to suppress the low-frequency noise originating from the 
surrounding nuclear spin flip-flops. This process, also 



termed nuclear spin diffusion, is the main source of de- 
coherence in paramagnetic systems containing magnetic 
nuclei, as is the present case of H@POSS. Therefore, the 
isolation of this decoherence mechanism can reveal the 
intrinsic spin-spin relaxation time which has to be com- 
pared with T 2 =160 [is, the longest coherence time of any 
solid state molecular radical that was recently reported 
for the 15 N@C6o system magnetically diluted in Ceo [3(. 
Moreover, by applying dynamical decoupling we identify 
the nuclear spin baths surrounding the unpaired electron 
spin and we probe decoherence processes that were pre- 
viously unexplored in paramagnetic trapped atoms. 

Although CPMG sequence is a well-established and 
frequently used method in NMR spectroscopy [14J, Il5| . 
its systematic application in EPR was only recently re- 
visited as a dynamical decoupling method in the con- 
text of spin-based quantum computing [ll|, [l6| ■ This is 
partly because paramagnetic systems have often broad 
EPR spectra that cannot be fully excited by microwave 



pulses leading to peculiarities 17[. Atomic hydrogen is 
an ideal testbed for applying CPMG decoupling because 
the solid state EPR spectrum exhibits narrow lines of the 
order of 0.1 mT. In order to ensure a small delocaliza- 
tion of electron spin density over the cage, we used the 
species with R=OSiMc3 shown in FigOJa), also known as 
QsMs [8]. Our previous study on the similar system with 
R=OSiMe2H revealed that for temperatures below 200 K 
the spin-spin relaxation is strongly enhanced by dynamic 
processes like rotation of the methyl groups [10( . There- 
fore, in our effort to avoid this additional decoherence 
process, we chose to perform decoupling experiments at 
T =200 K where the upper limit of Ti is virtually im- 
posed by the spin-lattice relaxation time T± =133 fis. 
FigHJb) shows that at this temperature the two-pulse 
echo decay is described by a stretch exponential function 



E(2t) = E exp 
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FIG. 1. (Color online) (a) Structure of the POSS molecule 
SigOi2(OSiMe3)8 with the encaged H atom, (b) Two-pulse 
electron spin echo decay measured at the low-field EPR tran- 
sition, Bo =318.6 mT, with superimposed stretched expo- 
nential fit using Eq. |l}. (c) CPMG decay measured at 
Bo =318.6 mT with two different r values and N = 30 pulses 
as a function of the time delay between the first 7r/2 pulse 
and the occurring echoes, t = 2r, At, ..., 2Nt. Inset: CPMG 
sequence used for dynamical decoupling. 
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FIG. 2. (Color online) Decay of the iV'th echo in a CPMG 
sequence using N = 10 as a function of interpulse delay t. 
(a) High-field EPR transition. Inset: pulse sequence showing 
the recorded echo, (b) Low-field EPR transition. The arrows 
in the enlarged time windows show the two r values used in 
the measurements of FigfJ] 



where r is the interpulse delay, n is a parameter deter- 
mined by the mechanism of phase memory decay and the 
rate, W, of the dephasing process relative to r, and Tm is 
the so-called phase memory time encompassing Ti and all 
other processes that cause electron spin dephasing 18]. 
At T =200 K we find T M =14.4 ^s and n =2.3, which im- 
plies a slow relaxation processes (Wt < 1) compatible 
with nuclear spin diffusion [l9j. 

FigOJc) demonstrates that the CPMG decay mea- 
sured with r = 3320 ns (blue squares) is accurately de- 
scribed by a monoexponential function with time con- 
stant Ti = 55.0 /^s. This is about four times longer 
than the value obtained with the two-pulse sequence, 
indicating a significant suppression of nuclear spin dif- 
fusion. Such dynamical decoupling is possible because 
the rate W of nuclear spin flips is of the order of their 
dipolar interaction (typically ~ 10 kHz for protons) cor- 
responding to a correlation time t c ^16 /is which is much 
larger than the used r value. Therefore, for short inter- 
pulse delays the decoherence process behaves like time- 
independent. Based on this concept, it is natural to as- 



sume that the longest relaxation decay corresponds to 
the minimum possible r value whose lower bound is equal 
to the spectrometer deadtime, Td- However, the CPMG 
decay obtained with r = 536 ns (FigJT] (c) red circles) 
shows clearly two peculiarities: first, the initial decay 
is surprisingly fast, even faster than the two-pulse echo 
decay; second, the total decay is not monotonic and can- 
not be described by a general stretch exponential. To 
the best of our knowledge, the dependence of the CPMG 
time constant on the time between the refocusing pulses 
has been documented once for an organic radical diluted 
in a mixture of deuterated solvents [13] , and has been 
correlated to the electron spin echo envelope modulation 
(ESEEM) effect ^. Moreover, by using CPMG pulse 
sequences on the nitrogen-vacancy center an anomalous 
decoherence effect has been recently reported 21 1. 

The above unusual behavior can be explained if we 
consider both the filter function of the applied CPMG 
pulse sequence, /(£; lj) = J f(t, i')e~ lwt dt', and the noise 



spectrum S(tu) of the system under study [22J. The de- 
coherence rate is determined by the overlap of these two 
functions. For a minimum overlap the coherence is pre- 
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FIG. 3. (Color online) (a) Baseline-corrected CPMG decays obtained with N =10 (blue traces) and their numerical simulations 
(red traces) for the two EPR transitions, (b) Absolute-value spectra obtained after FT of the time-domain experimental (blue 
lines) and simulated (red lines) signals. Asterisks mark "folded-back" peaks corresponding to even harmonics of the proton 
nuclear Zeeman frequency uh that are larger than the Nyquist frequency, (c) Noise spectral densities occurred from numerical 
simulations. Dashed lines denote the 1 H Larmor frequency. 



served for a long time, whereas for a maximum overlap 
the fastest decay occurs. Compared to the two-pulse se- 
quence the filter function of the CPMG exhibits peaks 
at higher frequencies. Therefore, apart from suppress- 
ing low-frequency processes, it can also probe an existing 
high-frequency noise that is otherwise 'silent' in a Hahn 
echo experiment. 

In order to understand the details of this peculiarity 
we performed CPMG experiments by measuring the in- 
tensity of the last echo of the pulse train as a function 
of t. This can give further insight into the spin dynam- 
ics of our system because each r value corresponds to 
a different filter function of the CPMG sequence. Fig|2] 
demonstrates that in addition to the exponential decay 
the echo shows sharp dips with a modulation depth that 
exceeds 0.5 for some r values. Applying dynamical de- 
coupling with these r values (e.g. r=536 ns as indicated 
by the up arrow of Fig(2fb)) result in very fast CPMG 
decays (see red circles in Fig FT]) whereas measurements at 
maximum echo positions give exponential CPMG decays 
with time constants close to 60 /is. 

The origin of the deep echo modulations can be best 
explored by analyzing the Fourier transform (FT) spectra 
shown in FigU^b) which contain peaks close to even har- 



monics of the 1 H and 29 Si nuclear Zeeman frequencies. 
These modulations are typical in ESEEM spectroscopy 
and are attributed to weak anisotropic hyperfine cou- 
plings between the unpaired electron and matrix mag- 
netic nuclei. Indeed, previous EPR studies have revealed 
such interactions both with protons of the cage organic 
substituents 12311 as well as with 29 Si nuclei of the core 
POSS cage [Irj. However, these couplings are weak and 
lead to shallow modulation depths as can be justified 
from the absence of dips in the two-pulse echo decay. 
Interestingly, the large number N of refocusing pulses in 
the CPMG sequence enhances the modulation depth con- 
siderably. Similar behaviors have been recently observed 
on nitrogen- vacancy centers (5 = 1) [24426J where the 
CPMG sequence was used to detect very weak hyper- 
fine couplings from remote 13 C nuclear spins. In that 
case, and for large external magnetic fields (weak cou- 
pling regime, wi 3> -B), the basic modulation frequency 
is given by uj± = \2u>i ± A\, where uj is the nuclear 
Zeeman frequency and A, B are the secular and pseu- 
dosecular parts of the hyperfine interaction, respectively. 
This frequency is equal to the sum of luq = \uji\ and 
u ±1 = i/(w/ ± A) 2 + B 2 , the corresponding ESEEM fre- 
quencies for 777,5 — and 7715 = ±1, in the limit of very 



small B [20(. For S = 1/2, which is the case of our 
present study, the sum of the two ESEEM frequencies 
uj a , ujp is given by [20f 



c 1.0 



Si 0.8 



uj + = 2|w/| 



\jQi\B 2 
Auj - A 2 ' 



(2) 



Contrary to S = 1, the weak coupling regime for S = 1/2 
results in sum frequencies that are virtually free from hy- 
pcrfinc coupling constants. This is the case for the weakly 
coupled protons of the cage organic groups for which the 
second term of Eq. ([2]) is estimated to be less than 5 kHz. 
For these protons the modulation frequencies match ex- 
actly 2vh and its higher harmonics, as can be seen from 
FigEJb). On the other hand, a careful observation at the 
low-frequency part of spectra reveals that the first silicon 
peak is shifted by about 75 kHz from 2v$i whereas for the 
higher harmonics the shifts are - as expected - more pro- 
nounced. We therefore deduce that the silicon hypcrfinc 
coupling is comparable to Vsi , in line with the previously 
obtained parameters \A\ =2.06 MHz and |B| =0.45 MHz 
for the inner silicon atom of the cage [10| . 

In order to quantify these observations we performed 
numerical calculations using two different approaches: 
for 29 Si we followed a complete simulation of the pulse 
sequence, whereas for the *H nuclear spin bath the cal- 
culation was done considering the filter function of the 
CPMG sequence and a Gaussian function for the noise 
spectral density [22j. FigJ3][a) and (b) show that both 
time and frequency domain data are in very good agree- 
ment with the simulated ones. For protons, where the 
free fitting parameters were the Gaussian function con- 
stants, the inferred noise spectral densities are centered 
at the X H Larmor frequency, have a maximum of about 
0.46 AiT/vTTi and a broadening of 15 kHz (FigJSfc)). 
From these Gaussian functions we calculate an rms mag- 
netic field noise 



'fS B (w)du =0.09 mT which is in 
good agreement with the observed EPR linewidth of 0.1 
mT. 

The validity of our model was also tested by simulat- 
ing CPMG decays for different r values. Fig|4]shows that 
the simulations reproduce well the experimental data and 
this further supports the argument that it is the proton 
high-frequency noise that accounts for the fast decoher- 
ence under special conditions of dynamical decoupling. 
On the other hand, a proper choice of r suppresses both 
low- and high-frequency noise leading to the longest time 
constant of about 56 /us (FiglU squares). This value is 
smaller than T\ and the expected Ti for the electron spin 
concentration of the sample under study [22[. Moreover, 
it is virtually constant for proper r values fulfilling 400 
ns< r < 4000 ns, indicating a white-noise behavior for 
frequencies between 60 kHz and 600 kHz, similar to a 
"frozen-core" effect predicted for Si:P 



27] 



Future sub- 
stitution of magnetic nuclei of the POSS cage (e.g. by 
perdeuteration) will elucidate whether this noise belongs 
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FIG. 4. (Color online) CPMG decays measured with differ- 
ent r values and N = 30 pulses. Open symbols, experimental 
data; filled symbols, simulations, (a) High-field EPR transi- 
tion, Bo =369.5 mT. (b) Low-filed EPR transition, B =318.6 
mT. Black lines correspond to monoexponential fits with time 
constants T2. 



to the tail of the proton low-frequency spectral density. 
In conclusion, the application of the CPMG sequence 
on H@QsM8 has shown that besides nuclear spin diffu- 
sion there are also two other sources of decoherence re- 
lated to the nuclear spin bath: first, the weakly coupled 
protons of the organic groups of the cage behave like a 
fluctuating magnetic field with the proton Larmor fre- 
quency, vh- This field induces a random phase on the 
central electron spin and creates additional dephasing. 
Second, for molecules possessing one or more 29 Si nu- 
clei in the core of the cage, the electron spin coherence 
is additionally modulated because the anisotropic hyper- 
fine interaction is of the order of the nuclear Zeeman 
interaction and thus induces a degree of entanglement 
between the electron and the 29 Si nuclear spin. These 
two sources of decoherence can be described as two dif- 
ferent limits of the electron-nuclear hyperfine coupling 
2J, |28[ but are not manifested in a typical Hahn echo 
experiment. Our work shows that although dynamical 
decoupling with multi-pulse sequences can significantly 
suppress low-frequency noise, it may also enhance deco- 
herence if the system bears a source of high-frequency 
noise; this should have high impact on quantum comput- 
ing schemes based on dynamical decoupling and solid- 
state spin qubits. Moreover, the increase of the phase 
memory time up to Tm=56 fj,s in a chemically unmodi- 
fied environment, poses atomic hydrogen encapsulated in 
POSS cages as a promising hybrid electron-nuclear spin 
system for quantum computing applications, comparable 
to endohedral fullerenes. 
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1 Theory and Simulations 

The numerical calculation of the CPMG signal can in principle be done using the spin Hamiltonian of the 

electron-nuclear coupled spin system under the unitary transformations of the pulse train. While this is 

^D ■ straightforward for one 29 Si nucleus, it becomes virtually impossible for 1 H nuclear spins due to their large 

■^ number and the wide distribution of their magnetic parameters. Therefore, for the latter case we assumed 

qJ a Gaussian function for the noise spectral density of the *H nuclear spin bath and the calculation of the 

<^ ■ CPMG decay was done by considering the filter function of the pulse train. 

1.1 Echo modulations due to coupling with a single inner-cage 29 Si spin 
The rotating frame spin Hamiltonian is given by [l| 

u = n s s z + uih + as z i z + bs z i x , (i) 



m 



43 



CZ2 

where O5 = ojs — oj mw is the offset of the electron Zeeman frequency ujs = g(3 e Bo/h from the mw frequency 
ujmw, wi = —g n (3 n Bo/h is the nuclear Zeeman frequency, g and g n are the electron and nuclear g-factors, 

C3 ! Pe and Pn are the Bohr and 29 Si nuclear magnetons, Bq is the static magnetic field along z-axis, and A, B 
describe the secular and pseudo-secular part of the hyperfine coupling. For an axially symmetric hyperfine 
_L ■ coupling tensor, L4_l,^4_l, ^.||]=[oiso — T, ai so — T, Oi so + 2T], where aj so and T are the isotropic the anisotropic 

£^ \ part of the hyperfine interaction, we have 

o 

U, ; A = a iso + T(3cos 2 6» - 1), B = 3Tcos6>sin6>, (2) 

with 6 being the angle between the static magnetic field Bq and the z— principal direction of the hyperfine 
^T ' tensor. Starting from the equilibrium density operator a eq = —S z , the first (n/2) x pulse creates electron 
,__! ■ spin coherence <To = R^ 1 (ft /2)a oq R x (Tr /2) , where R x (tt/2) = exp[—iir/2S x ]- After the iV'th {ft) y — pulse the 

echo intensity is given by 

{o- x )si = Re Tv(S x a N ), (3) 

^ ; 

^^ , where o~n is the density matrix at the time of the N th echo formation, t = 2Nt, 

CO 

a N = U^R-^ftp^aN-xUrRy^Ur, (4) 

, and U T = exp[— H-Lqt] is the propagator of free evolution. 

From our previous HYSCORE study [2J we obtained ai so — —2 MHz and T ~ —0.6 MHz for the inner 
c3 29 Si cage atoms. We used these values as starting point for the simulations which included integration over 
all possible orientations. Refinement of the parameters gave a^ so = —2.2 MHz and T = —0.65 MHz that 
reproduced very well the position and intensity of peaks close to even harmonics of us% (see Fig. ISlj) . Note 
that the high-field spectrum is more appropriate for this comparison because most of the 29 Si peaks do not 
overlap with 1 H harmonics. 

1.2 Echo modulations due to coupling with the 1 H spin bath 

We treat the surrounding l H nuclear spin bath as a source of classical noise with spectral density S(uj), 
i.e. a fluctuating magnetic field whose frequency is peaked at the X H Larmor frequency vr and distributed 
according to S(oj). The coherence decay is given by [3J 

(a x ) H = exp(-x(t)) with x(t) = j =£ S(oj)\f(t; W )| 2 , (5) 
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FIG. SI: (Color online) Absolute value FT spectra of experimental (blue lines) and separately simulated 
(green- 29 Si, red- 1 !!) signals of the last echo of a CMPG sequence with TV" =10. (a) High-field EPR transition, 
Bq =369.5 mT. (b) Low-filed EPR transition, Bq =318.6 mT. Black dashed lines denote even harmonics of 
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FIG. S2: (Color online) CPMG sequence for N =4 and the corresponding function f(t,t'). 



where /(£; u>) is the filter function of the pulse sequence given by the Fourier transform of the function f(t, t') 
with respect to t', 

f(t;u)= I f(t,t')e-^'dt>. (6) 

for N =4. For 



/ f(t,t')e-^'dt'. 

For the CPMG sequence the function f(t,t') is the Heaviside step function shown in Fig. 
the 1 H spectral density we assumed a Gaussian distribution of width a, 



S{u) = A 2 



1 



o-V2tt 



exp 



(u - U L ) 
2a 2 



2l 



(7) 



where ujl is the 1 H Larmor frequency and A is the coupling strength parameter having units of angular 
frequency. Another useful representation of this noise is through the magnetic field sensitivity, Sb(oj), 
defined as 



S b (uj) = — 

7e 



2vr 



(8) 



where 7 e /27r = 28.025 x 10 9 Hz/T is the gyromagnetic ratio and Sb{u) has units of T/vHz. Fig. 3(c) in 
main text shows noise spectral densities in this representation. 

Fig. [S3] compares the filter functions of the two-pulse echo and CPMG sequence with A^ = 10 for the 
same time evolution t =4 /xs. Fig. IS37 a) shows that, compared to two-pulse, the CPMG sequence not only 
suppresses low- frequency noise but also exhibits peaks at the high- frequency region where the 1 H noise 
appears. When the filter and noise spectral density functions overlap (Fig. IS3T b)). additional decoherence 



occurs. 



1.3 Overall simulation of echo modulations 

The total echo modulation is calculated as the product (cr x )H{o'x)si taking into account the natural abun- 
dance 4.68% of 29 Si. Considering only the SisOi2 core of the POSS cage, the binomial distribution predicts 
that Co = 68.13% of the cages have no 29 Si atom. Therefore, these species contribute only to 1 H modula- 
tions. The rest c = 31.87% of the cages have one {c\ = 26.78%), two [pi = 4.61%), three (03 = 0.45%), four 
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FIG. S3: (Color online) (a) Comparison of the filter functions for t = 4/is between the two-pulse echo (blue 
thin line) and CPMG with N = 10 (red thick line), (b) Details of the filter functions in the high frequency 
region together with the simulated 1 H spectral density (orange dashed line) for the low-field EPR transition. 

(04 = 0.03%) etc. 29 Si atom(s) in the core and therefore contribute also to 29 Si modulations. Assuming that 
all eight 29 Si atoms are magnetically equivalent, the total echo modulation is given by 



(crx) =^2c j (a x ) H [{(Tx)siW 



(9) 



In our simulations we used only the first four terms of Eq. (up to j = 3) because the rest of them have 
negligible contribution. Fig. IS4I (green thick lines) shows the contribution of 29 Si modulations to the total 



i = 3024 ns 




2.5 
Interpulse delay, t/)js 

FIG. S4: (Color online) Baseline-corrected CPMG decays obtained with N =10 (blue traces) and their 
numerical simulations (red traces) for the two EPR transitions. Green traces are the contribution of 29 Si 
modulations to the total signal. Arrows mark the r values used in measurements shown in Fig. 4 

signal along with the experimental data. Setting A and a as free varying parameters, the noise spectral 
densities were obtained by minimizing the sum of the residuals. The results gave A = 1.53 ± 0.05 x 10 7 rad/s 
and o-/2vr = 14.8 ± 0.5 kHz for the high-field, and A = 1.45 ± 0.05 x 10 7 rad/s and a/2ir = 15.0 ± 0.5 kHz 
for the low- field EPR transitions. This broadening corresponds to static magnetic field Bq inhomogeneity 



of about 0.3 mT which is much larger than the instrumental magnetic field inhomogeneity (using the same 
amount of sample the resolution of the liquid-state EPR spectrum acquired with the same spectrometer was 
about 10 /iT). Therefore, we can assign this linewidth to inhomogeneous broadening of the proton nuclear 
spin spectrum. Moreover, the obtained a values are comparable to the dipolar couplings between methyl 
group protons, i.e. 20 kHz, showing that nuclear spin flip-flops driven by dipolar interactions are indeed 
quite probable. 

1.4 Discussion on the obtained coherence time T 2 = 56 /is 

Here we discuss the observed T 2 value as a possible result of other relaxation effects like instantaneous 
diffusion (ID) and electron spin-spin interactions. The estimated electron spin concentration of our sample 
is about C = 1.3 x 10 16 spins/cm 3 . Since the two EPR transitions are very well separated, the time constant 
due to ID is given by [l[ 

" C 77wS £ sm TT> ( 10 ) 



r 2( iD) 18^3 h 2 

where fio is the permeability of free space and # 2 the rotation angle of the refocusing pulse. Setting 62 = n 
and g = 2 gives T 2 (jrj) = 186 [is which is larger than both relaxation times T\ = 133 /xs and T 2 = 56 /is at 
this temperature (T = 200 K). Moreover, similar measurements on a sample of half concentration revealed 
the same T 2 . We can therefore safely conclude that ID is not the dominant relaxation mechanism that 
determines the CPMG decay. On the other hand, Ti-induced spectral diffusion is possible because due 
to the large hyperfine splitting we cannot excite both EPR transitions simultaneously. Although CPMG 
refocuses interactions between resonant (A) and non-resonant (B) electron spins, random flips of the B 
spins with relaxation time T\ will contribute to the echo decay. This contribution is both temperature and 
concentration dependent. Preliminary measurements at very low temperatures (T = 5 K) showed that the 
spin-lattice relaxation time becomes extremely long, T\ = 57 ± 2 s, however, the maximum observed CPMG 
time constant of T 2 ~ 60 fj,s remained practically unchanged. This is a strong indication that Ti-spectral 
diffusion effects do not determine this decay. 

Another mechanism that can contribute to decoherence is flip-flopping between neighboring electron spin 
pairs driven by dipole-dipole interactions. This includes transfer of magnetization between the central and 
a neighbor electron spin (spin diffusion) or change of the resonance frequency of the central electron spin 
due to flip-flops between two neighbors (T 2 -type spectral diffusion). For this to be significant, the dipolar 
coupling of the two spins involved has to be larger or at least of the same order of magnitude as the difference 
in their resonance frequencies. The dipole-dipole coupling between two electron spins is given (in frequency 
units) by [l( 

^0 9\9lfi e /-, o 2/j \ n -,\ 

dd = 4irh rf 2 ( " 12) ' ( n) 

where ri 2 is the distance between the two spins and 612 is the angle between the static magnetic field Bq 
and the inter-spin axis. Because the frequency difference of the two EPR transitions is very large (about 
1.42 GHz), the effective concentration of the spins that can undergo flip-flops is C/2 ~ 6.5 x 10 15 spins/cm 3 
that gives an average spin separation ri 2 = 65 nm. For this distance the average dipolar coupling frequency 
is about 90 Hz, much smaller than the inhomogeneous broadening of about 3 MHz of the EPR lines. We 
therefore conclude that flip-flopping mechanisms involving electron spins are highly suppressed in our system. 

2 Materials and methods 

2.1 Sample preparation 

Octakis (trimethylsilyloxy) silsesquioxane (SisOi 2 (OSiMe3)8 or QsMg, CAS 51777-38-9) was prepared fol- 
lowing a two-step process: 

Octaanion salt synthesis (1): The octaanion [Sis0 2 o 8_ ]-8[Me4N] + was prepared according to Hagiwara 
et al. 0]. The mixture of TEOS/TMAOH/H 2 0/EtOH in molar ratio 1/1/10/10 was left for vigorous 
stirring for 3 days. After this reaction time, the mixture concentrated by vacuum and cooled. Hydrated 
crystals of the octaanion were observed. These crystals were dehydrated after heating at 60 °C in high 
pressure vacuum, lg of the isolated product was dissolved in a mixture of pyridine/THF (150/100) which 
were pretreated in an ice bath. After 30 min stirring the unreacted products were removed under reduced 



pressure. A portion of HC1 was added to neutralize the pyridine excess. The precipitated salt was isolated 
with filtration becoming a clear solution. The liquid was concentrated and the octaanion salt resulted. 

[Si 8 02o][Si(CH 3 ) 3 ] 8 (2): Trimethylchlorosilane (22.0 ml, 0.189 mol) was added in hexane (120 ml) and 
the mixture was stirred for 30 min at °C. The flask equipped with a reflux condenser and an additional 
funnel aiming at adding the octaanion butanol solution dropwise (1.0 g octaanion in 50.0 ml butanol) under 
N2 during 30 minutes period. The reaction was left for stirring additional 30 minutes and after that the 
mixture was extracted with hexane/brine twice and the organic layer was recovered. The organic layer was 
dried over Na2S04. The hexane was concentrated to give a white solid in a good yield. Compound 2 was 
characterized by FT-IR and their characteristic peaks are: Si-O-Si at 1100 cm -1 , Si-CH3 at 1260 cm -1 [5J. 

Encapsulation of atomic hydrogen was performed with 7-irradiation in steps using a 60 Co source. The 
accumulated dose was measured with Red Perspex Dosimeters, Type 4034 AD. 54 mg of POSS powder were 
placed in a sealed vial filled with O2 (acting as radical scavenger) and irradiated for 20.9 days resulting in a 
total dose of 255 kGy. After 7-irradiation, 27.8 mg of the sample was transferred into an EPR quartz tube. 
Comparison of the cw EPR spectrum (in presence of O2, unsaturated conditions) with a standard sample 
gave an estimate of 1.3 x 10 16 spins/cm 3 for the electron spin concentration. To avoid paramagnetic oxygen 
in pulsed EPR measurements, the quartz tube was filled with He gas following the procedure described 
previously [2|. 



2.2 Pulsed EPR spectroscopy 

EPR measurements at X-band {oj mw /2'K = 9.698 GHz) were carried out on a Bruker ESP 380E spectrometer 
equipped with an EN 4118X-MD4 Bruker resonator. Measurements at T = 200 K were performed with a 
helium cryostat from Oxford Inc. The microwave frequency was measured with a HP 5350B microwave 
frequency counter and the temperature was stabilized with an Oxford ITC4 temperature controller. Two- 
pulse echo and GPMG measurements used t w / 2 = 16 ns and t v = 32 ns. For the CPMG sequence the 
maximum number of pulses N is limited by the spectrometer to 30. For CPMG experiments with constant r 
values the measurements were recorded with a HP Infinium 54810A oscilloscope which allowed the acquisition 
of the entire time trace in a single shot (see Fig. IS5p . All measurements were performed with repetition 
rates smaller than l/5Ti = 1.5 kHz in order to avoid saturation. 
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FIG. S5: (Color online) Oscilloscope-recorded CPMG time trace for N = 30 and r = 720 ns. Number of 
averaging measurements, 1024; total acquisition time, 30 s. The trace is the difference between on- and 
off-resonance signals in order to remove unwanted baseline artifacts induced by microwave pulses. 



2.3 Data manipulation 

The data were processed with the program MATLAB 7.0 (The MathWorks, Natick, MA). For r values 
corresponding to maximum echo intensities in CPMG measurements, T2 relaxation times were determined 
by fitting the time traces with single exponential functions, whereas for obtaining Tm from two-pulse echo 
decay a stretched exponential function was appropriate. The time traces of variable r CPMG experiments 
(Fig2. main text) were baseline corrected with a single exponential, apodized with a gaussian window and 
zero filled. After Fourier transform the absolute-value spectra were calculated. CPMG traces were simulated 
with a program written in-house based on the theory described above. 
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